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Abstract 

We examine the behavior of the variable Eddington factor for a relativisticahy moving radiative flow in 
the vertical direction. We adopt the "one-tau photo-oval" approximation in the comoving frame. Namely, 
the comoving observer sees radiation coming from a closed surface where the optical depth measured from 
the observer is unity; such a surface is called a one-tau photo- oval. In general, the radiative intensity 
emitted by the photo-oval is non-uniform and anisotropic. Furthermore, the photo-oval surface has a 
relative velocity with respect to the comoving observer, and therefore, the observed intensity suffers from 
the Doppler effect and aberration. In addition, the background intensity usually depends on the optical 
depth. All of these introduce the anisotropy to the radiation field observed by the comoving observer. As 
a result, the relativistic Eddington factor / generally depends on the optical depth r, the four velocity 
and the velocity gradient du/dr. In the case of a plane-parallel vertical flow, we found that the relativistic 
variable Eddington factor / generally decreases as the velocity gradient increases, but it increases as the 
velocity increases for some case. When the comoving radiation field is uniform, it is well approximated by 
3/ ~ 1/[1 + (16/15)(— (iix/7(ir) + (— (ii^/7(ir)^'^~^]. When the radiation field in the inertial frame is uniform, 
on the other hand, it is expressed as / = (1 + 3;^^)/ (3 + These relativistic variable Eddington factors 
can be used in various relativistic radiatively-driven flows, such as black- hole accretion flows, relativistic 
astrophysical jets and outflows, and relativistic explosions like gamma-ray bursts. 

Key words: accretion, accretion disks — astrophysical jets — gamma-ray bursts — radiative transfer 
— relativity 



1. Relativistic Variable Eddington Factor 

In the moment formalism of (relativistic) radiation hy- 
drodynamics, in order to close the moment equations 
truncated at a finite order, we need a closure rela- 
tion^ such as the Eddington approximation or more com- 
plex expressions. In non-relativistic static atmospheres 
(Chandrasekhar 1960; Mihalas 1970; Rybicki & Lightman 
1979; Mihalas & Mihalas 1984; Shu 1991; Peraiah 2002; 
Castor 2004), for example, many researchers have used the 
Eddington approximation and variable Eddington factors 
as closure relations (Milne 1921; Eddington 1926; Kosirev 
1934; Chandrasekhar 1934; Hummer & Rybicki 1971; 
Wilson et al. 1972; Tamazawa et al. 1975; Unno & Kondo 
1976; Masaki & Unno 1978). As is well known, the stan- 
dard Eddington approximation with the Eddington factor 
of 1/3 is valid when the radiation field is almost isotropic. 
In the spherically symmetric case, where the radiation 
field becomes anisotropic towards the outer optically thin 
region, they often used the variable Eddington factor, 
which depends on the optical depth (e.g., Tamazawa et 
al. 1975). 

In relativisticahy moving flows (Mihalas et al. 1975 
Mihalas et al. 1976a, b; Thorne 1981; Thorne et al. 1981 
Flammang 1982, 1984; Mihalas & Mihalas 1984; Nobih 
et al. 1991, 1993; Kato et al. 1998, 2008; Castor 2004; 
Mihalas & Auer 2001; Park 2001, 2006; Takahashi 2007), 



on the other hand, the standard Eddington approxima- 
tion is adopted in the comoving frame, and then trans- 
formed to the inertial frame, if necessary (Castor 1972; 
Hsieh & Spiegel 1976; Fukue et al. 1985; Sen & Wilson 
1993; Baschek et al. 1995, 1997; Kato et al. 1998, 2008). 
Even for a comoving observer, however, the comoving ra- 
diation field may become anisotropic due to various rea- 
sons as shown below. Instead of the standard Eddington 
approximation in the comoving frame, for such a rela- 
tivistic case several types of a variable Eddington factor^ 
which depends on the velocity and its gradient as well as 
the optical depth, were proposed (Fukue 2006; Fukue & 
Akizuki 2006, 2007; Akizuki & Fukue 2008; Fukue 2008a, 
b; Koizumi & Umemura 2008). However, the moment 
formalism of relativistic radiation hydrodynamics is yet 
incomplete^ in a sense that we do not have the sufficiently 
adequate closure relation. 

There are mainly three reasons, which cause the 
anisotropy of the comoving radiation field in the relativis- 
ticahy moving flow. (1) The optical depth effect. In the 
outer region where the flow becomes optically thin, sim- 
ilar to the non-relativistic case, the comoving radiation 
field would be anisotropic. (2) The velocity gradient ef- 
fect. When the radiative flow is accelerated up to the 
relativistic regime and there is a strong velocity gradi- 
ent in the direction of the flow, the velocity fields as well 
as the density distribution are no longer uniform even in 
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the comoving frame. Then, the comoving radiation field 
also becomes non-uniform. Furthermore, due to the rel- 
ative speed between the observer and the radiation field, 
the observed radiation field suffers from the Doppler ef- 
fect. As a result, the comoving radiation field becomes 
anisotropic. (3) The effect of the relativistic speed itself. 
Although we do not know the radiative intensity in the 
inert ial frame, in the highly relativistic regime it is gen- 
erally affected by the strong aberration effect when it is 
converted to the comoving frame. As a result, except for 
some special cases, the comoving radiative intensity gen- 
erally has strong anisotropy. 

Hence, the relativistic variable Eddington factor 
(RVEF) generally depends on the optical depth, the fiow 
velocity, and the velocity gradient. In Fukue (2006) and 
Akizuki and Fukue (2008), they proposed such variable 
Eddington factors from physical view points in the plane- 
parallel and spherical cases, respectively. In Fukue (2008a, 
b), on the other hand, he derived semi-analytically vari- 
able Eddington factors for the plane-parallel fiows in the 
vertical direction, and found that the Eddington factor de- 
creases as the velocity-gradient increases. In these studies 
(Fukue 2008a, b), however, the treatments were rather re- 
strictive; e.g., the comoving radiation field is assumed to 
be uniform. Thus, in the present paper we extend the 
previous work under the similar procedure to examine the 
above three reasons more extensively. 

In the next section we describe the shape of the sur- 
face (one-tau photo-oval), where the optical depth mea- 
sured by the comoving observer is unity, and derive the 
expression of the one-tau photo-oval. In section 3 we nu- 
merically calculate the comoving radiation field within the 
photo-oval, and the relativistic variable Eddington factor 
for several cases. The final section is devoted to conclud- 
ing remarks. 

2. One-Tau Photo-Oval and Photo- Vessel 

Let us suppose a relativistic radiative fiow, which is 
accelerated in the vertical (z) direction, and a comoving 
observer, who moves upward with the fiow (figure 1). 

Sufficiently deep inside the flow, where the optical depth 
is very large, the mean free path of photons is very short. 
In such a case, within the range of the mean free path, 
for the comoving observer the flow is seen to be almost 
uniform; the velocity gradient and resultant density gra- 
dient are negligible. Hence, the mean free path is the 
same in all directions, and the shape of the surface where 
the optical depth measured from the comoving observer is 
unity, is almost a sphere; we call it a one-tau photo-sphere 
(a dashed circle in flgure 1). As a result, the comoving 
radiation fleld is nearly isotropic, and a usual Eddington 
approximation in the comoving frame is valid. 

On the other hand, if the velocity gradient becomes very 
large, the behavior of the radiation fleld may be affected 
by changes in the hydrodynamics of the flow even on a 
lengthscale comparable to that of the photon mean free 
path. Hence, the mean free path becomes longer in the 
downstream direction than in the upstream and other di- 



surface 



Z 



b^se 



photo- 
vessel 



photo- 
oval 



photo- 
sphere 



Fig. 1. Schematic picture of a relativistic radiative flow in 
the vertical direction. The flow is accelerated in the vertical 
(z) direction, and has a velocity gradient. The dashed curves 
are one-tau photo-ovals observed by a comoving observer. 

rections, and the shape of the one-tau range elongates in 
the downstream direction; we call it a one-tau photo-oval 
(a dashed oval in flgure 1). As a result, the comoving ra- 
diation fleld becomes anisotropic due to various reasons 
as described in the introduction, and we should modify 
the usual Eddington approximation. 

In addition, when the optical depth is sufliciently small 
and/or the velocity gradient is sufliently large, the mean 
free path of photons in the downstream direction become 
less than unity. Hence, the shape of the one-tau range 
is open in the downstream direction; we call it a one-tau 
photo-vessel (a dashed hemi-circle in flgure 1). 

In order to obtain an appropriate form of the relativistic 
variable Eddington factor in these regimes, we thus care- 
fully treat and examine the radiation fleld in the comoving 
frame. 

In previous papers (Fukue 2008a, b), we have examined 
the one-tau photo-oval, derived the variable Eddington 
factor semi-analytically, and proved that the Eddington 
factor decreases with the velocity gradient. In these stud- 
ies, however, the treatment were rather restricted, and in 
the present study, we thus examine the behavior of the rel- 
ativistic variable Eddington factor for more general cases, 
and seek the appropriate form of the Eddington factor 
f{r^u^du/dr)^ where r is the optical depth and u (=7/3) 
is the four velocity, /3 and 7 (= l/y^l — P'^) being the nor- 
malized velocity and the Lorentz factor, respectively. 

We flrst derive the shape of the one-tau photo-oval using 
the optical depth r, while in the previous papers we have 
used the vertical coordinate z. 

We deflne the one-tau photo-oval as the surface where 
the optical depth measured by the comoving observer is 
unity. The situation is schematically illustrated in flgure 
2. We assume that the comoving observer in the vertical 
flow is located at z = or r = tq, where the flow four 
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Fig. 2. One-tau photo-oval around a comoving observer in 
the vertical (z) one-dimensional radiative flow. The comoving 
observer is located at r = tq, where the flow four speed is 
u = uq. In the s-direction, that forms an angle with the 
downstream direction, the mean free path is set to be ^o- 

speed is u = Uq. In the s-direction, that forms an angle 
with the downstream direction, the mean free path of 
photons is /q- The relation among these quantities is 

z — Zq = scosO. (1) 

The continuity equation for the stationary, one- 
dimensional relativistic flow is 

pcu = J (= const.), (2) 

where p is the proper gas density, u the four velocity, and 
J the mass-flow rate per unit area. In addition, the optical 
depth r is defined by 



dr = —npdz = —npdscosO^ 



(3) 



where is the opacity, which is assumed to be constant 
in the present analysis, and we used equation (1). 

In this paper we use a linear approximation for the flow 
field; that is to say, around the position of the comoving 
observer the flow four speed is expanded as 

du , . 

u = uo^-r- [r-To), (4) 
dr 

where du/dr\o is assumed to be constant. 

Using equation (3), the optical depth Tg along the s- 
direction is expressed as 

Ts = I npds = - I ——^ — TTTT- (^) 



cos 6 cos 6 

Thus, the shape of the one-tau photo-oval, where = 1, 
is finally determined by the condition: 



Obviously, for tq > 1 the one-tau region is close to be 
a photo-oval, whereas it is open towards a downstream 
direction to be a photo- vessel for tq < 1. 

3. Comoving Radiation Field and the Relativistic 
Variable Eddington Factor 

We can now compute the radiation field received by 
the comoving observer, and derive an expression for the 
Eddington factor in the comoving frame (figure 3). 

In a static and optically thick atmosphere, the radia- 
tion field is isotropic and uniform. In the present moving 
atmosphere, on the other hand, there are three reasons 
for which the radiation field is not uniform and isotropic, 
as described in the introduction; the optical depth effect, 
the velocity gradient effect, and the effect of the relativis- 
tic speed itself. 

First, the radiative intensity / in the comoving frame 
emitted from the one-tau photo-oval walls is generally a 
function of r and p. In general the radiative intensity 
increases with the optical depth, and therefore, the in- 
tensity from the upstream direction is slightly larger than 
that from the downstream direction in the present one- 
dimensional vertical flow. This effect of non-uniformity of 
intensity generally acts as a force to accelerate a comoving 
observer. This optical depth effect is prominent for small 
optical depth, similar to the static atmosphere. 

Second, although in the usual static atmosphere the p- 
dependence is safely ignored for large optical depth, this is 
not the present relativistic case. Namely, in the relativis- 
tic flow there appear the aberration and Doppler effects 
between the inert ial and comoving frames. Hence, except 
for very special cases, the comoving radiative intensity 
would have strong anisotropy. 

Thirdly, if there is a velocity gradient in the flow, the 
intensity observed by the comoving observer is redshifted 
(Doppler shifted) due to the velocity difference between 
the comoving observer and the one-tau photo-oval walls. 
In an accelerating flow, where the flow speed increases 
toward the downstream direction, the relative velocity is 
generally positive (figure 3) except for some special direc- 
tion {6 = 7r/2). Hence, the Doppler shift of intensity also 
causes anisotropy of the radiation field at the position of 
the comoving observer. 

Neglecting the local emissivity, the frequency-integrated 
intensity /co observed in the comoving frame at the 
position To of the comoving observer is related to the 
frequency-integrated intensity /(r, p) emitted at the one- 
tau photo-oval and to the redshift z by 



(i + z)4 



(7) 



The relative speed A/3 between the comoving observer and 
the one-tau photo-oval walls is given by the relativistic 
summation law as 



A/3-- 



13 -Po 



(8) 



r = To — cos^. 



(6) 

where 
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Fig. 3. Radiation field in the one-tau photo-oval around a 
comoving observer in the vertical one-dimensional radiative 
fiow. The radiation field may become anisotropic, since in 
general the emitted intensity is not uniform and isotropic, 
it is redshifted due to the velocity difference, and there is 
aberration between the comoving and inert ial frames. 



Uo 



(9) 
(10) 



It should be noted that as long as the velocity gradient is 
positive (accelerating flow) the relative speed (8) is always 
positive (or zero), and the one-tau photo-oval is seen to 
expand (figure 3). 

Using the relative speed above, the redshift z is ex- 
pressed as 

1 + A;^cos6> 



1 



(11) 



Now, all the quantities are calculated for given tq, uq^ 
du/dr\o^ and I{r^fi). Once the observed intensity /co is ob- 
tained by equation (7), the radiation energy density £^co, 
the radiative flux Fco, and the radiation pressure Pco niea- 
sured by the comoving observer are calculated respectively 



cEc. 



I 



(12) 



CPno = / ho COS^ OdVtn 



(14) 



Finally, the relativistic variable Eddinton factor in the 
comoving frame is given by 

f^^- (15) 



In the following subsections, we consider in turn several 
cases; the velocity gradient effect, the optical depth effect, 
and the strong anisotropy of the radiative intensity. 

3.1. Uniform Case with Photo-Vessel 

We first consider the velocity gradient effect, similar 
to the previous studies (Fukue 2008a, b), and also the 
photo- vessel case, which was not considered in the previ- 
ous studies. To clarify the velocity gradient effect, we here 
assume that the comoving radiation field is uniform: 



/(r,/i) = / = const. 



(16) 



The numerical results of this uniform case are shown 
in figure 4 for small uq and figure 5 for large ixo- In fig- 
ures 4a and 5a the relativistic variable Eddington factors 
f{u^ —du/dr) multiplied by 3 are plotted in the parameter 
space, whereas they are plotted for several fixed uq with 
fitting curves in figures 4b and 5b. 

The behavior of RVEF for small uq is shown in figure 
4 (cf. Fukue 2008a). In this case RVEF decreases, as the 
velocity gradient becomes large, as analytically proved in 
Fukue (2008a). In figure 4b for fixed uq (= 0.01,0.05,0.1) 
the values of 3/ are plotted as a function of — du/dr \o 
by solid curves; all curves overlap each other. One fit- 
ting curve is also plotted by a thick dashed one, and the 
numerical results are well fitted by 



3/' 



. 16 

IH V- 

15^ 



where 

y = -du/dT\o, 



(17) 



(18) 



within an error of a few percent. It should be noted that 
we chosed the coefficient, 16/15, so that theis fitting equa- 
tion reduces to 

3/ = i-^y (19) 

in the linear regime as analytically proved in Fukue 
(2008a). 

The behavior of RVEF for large uq is shown in figure 5 
(cf. Fukue 2008b). In this case RVEF decreases, as the 
logarithmic velocity gradient becomes large, as shown in 
Fukue (2008b). In figure 5b for fixed uq (= 1,5,10) the 
values of 3/ are plotted as a function of —dlnu/dr\o by 
solid curves; except for the marginally case of txq = 1 other 
curves overlap each other. One fitting curve is also plotted 
by a thick dashed one, and the numerical results are well 
fitted by 



(13) 3/^1 



16 
15' 



:a-\- a 



1.5 



where 



-d\nu/dr\o. 



(20) 



(21) 



Although in Fukue (2008b) we have used an exponential 
fitting curve (the dotted curve in figure 5b), the present 
power law type is well fitted at large a. 
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Fig. 4. Behavior of the relativistic variable Eddington factor 
f(u, —du/dr) multiphed by 3 in the case of small u. (a) 3/ are 
plotted as a function of the four velocity no {x axis) and the 
four- velocity gradient —(in/dr |o (y axis), (b) 3/ are plotted 
as a function of the four- velocity gradient —du/dT\Q for four 
fixed UQ {= 0.01,0.05,0.1; solid curves), and one fitting curve 
is also shown (dashed curve). 



Fig. 5. Behavior of the relativistic variable Eddington fac- 
tor /(n, —du/dr) multiplied by 3 in the case of large u. (a) 3/ 
are plotted as a function of the four velocity uq {x axis) and 
the logarithmic four- velocity gradient —d\nu/dT\Q (a axis), 
(b) 3/ are plotted as a function of the logarithmic four-ve- 
locity gradient —d\nu/dT\o for three fixed uq {= 1,5,10; solid 
curves) from top to bottom, and one fitting curve is also shown 
(dashed one), with an exponential curve (dotted curve). 
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Mixing the small and large uq cases, and considering 
several other patterns, the relativistic variable Eddington 
factor in the uniform field is well fitted by 



3/' 



^+15 



1 du 
7 dr 



ldu\ 
J dr J 



1.6-2 



(22) 



We here assumed that the comoving intensity is uni- 
form, the above result is only due to the redshift effect, 
similar to the previous studies (Fukue 2008a, b). 

Here, we briefly examine the case of small optical depth, 
a photo- vessel (figure 1). That is, as was stated, the one- 
tau region is open towards a downstream direction to be 
a photo- vessel for tq < 1 [equation (6)]. 

The behavior of RVEF for small r is shown in figure 
6. In figure 6a the values of 3/ are plotted as a func- 
tion of —du/drlo for several values of tq by dashed curves 
(tq = 0.1,0.2,0.3,0.4,0.5 from bottom to top) and by solid 
curves (tq = 0.6,0.7,0.8,0.9,1.0 from bottom to top) for 
U() = 0.1. The tendencies are similar, but there exist weak 
deviations. 

In figure 6b the value of 3/ is plotted as a function of 
To for 1^0 = 0.1 and —du/dr\o = 0. Even if there is no 
velocity gradient, the value of RVEF changes with the 
optical depth. This is understood as follows. When the 
optical depth tq is larger than unity (photo-oval) , / = 1 /3 
in the uniform radiation field without the velocity gradi- 
ent. When the optical depth becomes smaller than unity 
(photo- vessel), the Eddington factor also becomes smaller 
than 1/3, since the radiation field is no longer isotropic; 
the radiation from the downstream direction vanishes. 
When the optical depth approaches zero, however, the 
Eddington factor recovers to 1/3, since the Eddington fac- 
tor above the infinite plane is 1/3. The dependence on the 
optical depth in figure 6b is well fitted by 



1 



1 



(23) 



which is shown by a thick dashed curve in figure 6b, with- 
out any significant error. 

3.2. Weakly Non-Uniform Case 

Next, we examine the optical depth effect, where the 
background radiation field is not uniform. We assume 
that the comoving radiation field is the Milne-Eddington 
type: 



3 _ / 2 

I{T,ll) = -liT+- 



(24) 



Here, r is the optical depth, which is related to the ob- 
server's quantities by equation (6), and fi is the direction 
cosine in the comoving frame, which is related to the ob- 
server's quantities by 



(25) 



if we consider the aberration effect between the observer 
and the radiation site on the photo-oval in the comoving 
frame. 



5/10'; 
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Fig. 6. Behavior of the relativistic variable Eddington fac- 
tor f(u,—du/dr) multiplied by 3 for small r. (a) The depen- 
dence on the velocity gradient is shown for several values of tq 
(= 0.1,0.2,0.3,0.4,0.5 from bottom to top of dashed curves) 
(=0.6,0.7,0.8,0.9,1.0 from bottom to top of solid curves) for 
wo = 0.1. (b) The dependence on the optical depth is shown 
in the case of uq = 0.1 and —du/dr\o = 0. One fitting curve 
is also shown (dashed curve). 
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3.3. Strongly Anisotropic Case 
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Fig. 7. Behavior of the relativistic variable Eddington fac- 
tor f(u,—du/dr) multiplied by 3. It is plotted as a function 
of the four velocity uq {x axis) and the four-velocity gradi- 
ent —du/dr\o (y axis) in (a) the uniform case and in (b) the 
Milne-Eddington case for tq = 1. 

Using this intensity distribution in the comoving frame, 
we easily compute the radiation field and the Eddington 
factor. The numerical results are shown in figure 7. In fig- 
ures 7a and 7b the relativistic variable Eddington factors 
f{u^—du/dr) multiplied by 3 are plotted in the parameter 
space for the uniform case (figure 7a) and for the Milne- 
Eddington case (figure 7b). 

As is seen in figure 7, the distribution of RVEF of 
the Milne-Eddington case is not so much different from 
that of the uniform case. For larger uq the difference be- 
comes much small. Hence, the weak non-uniformity or 
weak anisotropy does not significantly affect the values of 
RVEF. 



Thirdly, we investigate the effect of the relativistic speed 
itself. Although there may be various kinds of anisotropy, 
we consider the simple limiting case; we assume that the 
radiation field in the inertial frame is uniform. Then, the 
comoving radiation field is expressed as 

/(r,/i,/3) = ---^--^/, (26) 

where P is the flow speed in units of the speed of light, 
7 the Lorentz factor, and /i the direction cosine in the 
comoving frame given by equation (25). Thus, due to 
the aberration effect the comoving radiation field becomes 
strongly anistropic. Namely, when the radiative intensity 
in the comoving frame is uniform, the intensity in the 
inertial frame becomes anisotropic in the sense that the 
observer in the inertial frame receives the strong radiation 
from the upstream direction. When the radiative intensity 
in the inertial frame is uniform, on the other hand, the 
intensity in the comoving frame becomes anisotropic in 
the sense that the comoving observer receives the strong 
radiation from the downstream direction. 

The numerical results of this anisotropic case are shown 
in figure 8. In figure 8a the relativistic variable Eddington 
factor f{u^—du/dr) multiplied by 3 is plotted in the pa- 
rameter space, whereas it is plotted as a function of uq 
with fitting curves in figure 8b. 

The behavior of RVEF shown in figure 8a is quite im- 
pressive, since it seems that / does not depend on the ve- 
locity gradient, but depends only on the velocity. Indeed, 
in figure 8b the values of 3/ are plotted as a function of 
Uq for several fixed —du/dr\o (= 0,0.5,1) by solid curves, 
and all curves overlap each other. One fitting curve is also 
plotted by a thick dashed one, and the numerical results 
are well fitted by 



3/ = 



3(1 +41(2) 3(l + 3/?2) 



3 + 41(2 



3 + /32 



(27) 



As comparisons, RVEF proposed in Fukue (2006) are also 
plotted by the thin dashed one (3/ = 1 + 2/3^) and the thin 
dotted one (3/ = 1 + 2/3). 

As is seen in figure 8b, equation (27) reproduces the 
numerical results very well. This is interpreted as follows. 
When the radiation intensity in the inertial frame is uni- 
form, the radiation energy density ^lab, the radiative flux 
Flab, and the radiation pressure Piab in the inertial (labo- 
ratory) frame become, respectively. 



c£^lab ■ 



lab ■ 



:47r/, 
0, 

4-- 



(28) 
(29) 

(30) 



for To > 1. Hence, from the Lorentz transformation the 
radiation energy density Eco^ the radiative flux Fco, and 
the radiation pressure Pco in the comoving frame become, 
respectively, 

-2 (cEiab-2/3Fiab+/3'cPiab) 
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Fig. 8. Behavior of the relativistic variable Eddington factor 
f{u,—du/dr) multiphed by 3 in the strongly anisotropic case, 
(a) 3/ are plotted as a function of the four velocity uq (x axis) 
and the four- velocity gradient —du/dr\o {y axis), (b) 3/ are 
plotted as a function of the four velocity uq for three fixed 
du/drlo (= 0,0.5,1.0; solid curves), and one fitting curve is 
also shown (dashed curve), with two comparisons (thin dashed 
and dotted curves). 



Fco = 7' [(1 + /^')^lab - /^(C^lab + cPiab)] 



rP — 

CO — 



-4 n 

-47rJ-7'/3, 
72 {p^cE,,^, ■ 



2/3Fiab + cPiab) 



[Vol. , 
(31) 

(32) 

(33) 



Thus, the relativistic variable Eddington factor analyti- 
cally becomes 



Em 



1 



3 + /52 



(34) 



4. Concluding Remarks 



In this paper, we have derived the relativistic variable 
Eddington factor for a plane-parallel flow that is accel- 
erating in the vertical direction; we have introduced the 
one-tau photo-oval observed by the comoving observer, 
and then calculated the comoving radiation field and the 
Eddington factor for several situations, and found that 
the Eddington factor in the relativistic regime generally 
depend on the optical depth, the fiow velocity, and the 
velocity gradient. 

Of these, the dependence of RVEF on the optical depth 
is weak similar to the non-relativistic plane-parallel atmo- 
sphere. As for the dependence on the velocity gradient, 
on the other hand, RVEF generally decreases as the ve- 
locity gradient increases. When the velocity gradient is 
small, RVEF linearly decreases with the velocity gradient, 
as proved in Fukue (2008a). When the velocity gradient 
becomes large, RVEF decreases in the power law manner 
of the logarithmic velocity gradient (cf. Fukue 2008b). 
Finally, as for the dependence on the relativistic speed, 
RVEF can increases to be on the order of unity, if the 
comoving intensity is sufficiently anisotropic. 

We have considered the case of the plane-parallel verti- 
cal fiow. In the spherical fiow, on the other hand, there 
exists a geometrical dilution effect, and the situation can 
be somewhat different, as terms appear (e.g., Peraiah 
2002). The present approach, however, can be extended 
to treat also the spherical case and we will examine it in 
future work. 

The treatment of the relativistic Eddington factor pre- 
sented in this paper may turn out to be applicable in 
various aspects of relativistic astrophysics with radiation 
transfer; i.e., black- hole accretion ffows with supercriti- 
cal accretion rates, relativistic jets and winds driven by 
luminous central objects, relativistic explosions including 
gamma-ray bursts, neutrino transfers in supernova explo- 
sions, and various events occured in the proto universe 
(cf. Fukue 2008b). The treatment of the relativistic 
Eddington factor presented in this paper will help in inves- 
tigating these types of astrophysical problems that involve 
the solution of the fully relativistic radiation hydrodynam- 
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